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AF  FUNCTION  -  Weapon  delivery  and  defenses,  transport, 
advanced  air-breathing  engines. 


DEFICIENCY  -  Insufficient  under standing  of  the  basic  phys¬ 
ical,  chemical  and  fluid  dynamic  processes  of  ignition,  stabil¬ 
ity,  and  efficiency  of  catalytic  combustion.  Lack  of  guidelines 
for  predicting  catalytic  combustor  performance  and  for  solution 
of  existing  combustor  difficulties. 

OBJECTIVE  -  To  clarify  the  relative  importance,  and  to 
formulate  realistic  analytical  representation  of  homogeneous, 
heterogeneous  kinetics  and  transport  processes  in  catalytic 
combustion  phenomena  associated  with  advanced  air-breathing 
combustion  systems. 

HOW  WORK  CONTRIBUTES  -  Will  provide  additional  understand-  . 
ing  and  needed  realistic  analytical  modeling  of  homogeneous  and l^J  LrJ 
heterogeneous  high  temperature  catalytic  combustion  processes  ® 


not  now  available.  Will  contribute  to  establishing  realistic  •  l: 

guidelines  and  techniques  for  the  design  of  efficient,  stable 

jet  engine  catalytic  combustors.  gp 

APPROACH  -  Theoretical  and  experimental  studies  will  be 
made  of  basic  fluid  dynamic,  physical  and  chemical  processes  in 
catalytic  combustion  associated  with  air-breathing  propulsion 
systems.  The  relative  importance  of  gas  phase  kinetics,  heat 
transfer,  mass  diffusion,  and  surface  chemical  kinetics  will  be 
assessed.  The  practical  phenomena  will  be  experimentally  simu¬ 
lated.  Deficiencies  in  existing  mathematical  models  will  be 
demonstrated  and  improved  models  formulated  based  on  experimen¬ 
tal  data  and  field  observations.  Various  monolithic  and  packed- 
bed  catalyst  candidates  for  advanced  combustor  design  will  be 
studied  over  a  range  of  operating  conditions  characteristic  of 
advanced  air-breathing  propulsion  engines  with  various  hydro¬ 
carbon  fuels.  The  two-dimensional  model  for  laminar  and  turbu¬ 
lent  boundary  layers  with  multiple  gas  and  surface  reactions 
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pheric  pressure  over  platinum/alumina/cordierite  catalysts  were  carried  out 
-ver  a  broad  range  of  equivalence  ratios  and  flow  velocities.  Substrate  (honey¬ 
comb)  temperatures,  chemical  speciss  at  the  catalyst  were  measured.  A  substan¬ 
tial  data  base  for  catalytic  oxidation  C^H'g  over  palatinum  catalysts  was  obtain¬ 
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ed.  When  the  experimental  wadi  temperature  is  used  as  the  boundary  conditions. 
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STATUS  OF  CATALYTIC  COMBUSTION  STUDY 


During  the  period  covered  by  this  report,  lean  combustion 
of  carbon  monoxide  in  air  and  in  inert/oxygen  mixtures,  and  of 
propane  in  air  was  modeled,  and  the  predictions  were  compared 
with  the  experimental  data  already  obtained.  A  new  catalyst 
material,  integral  with  the  substrate,  was  also  designed  and 
is  in  the  process  of  being  tested. 

A.  EXPERIMENTAL  APPARATUS 


The  new  experimental  apparatus  used  to  collect  the  data  is 
shown  in  Fig.  1.  With  respect  to  the  old  test  rig,  described 
in  the  previous  (6/1/79  -  5/31/80)  Annual  Technical  Report,  the 
following  modifications  have  been  made:  a)  an  Ingersol 1-Rand 
compressor  has  been  installed  that  can  be  connected  to  the  air 
feedline,  substituting  for  the  Atlas-Copco  compressor.  The 
maximum  pressure  can  now  be  increased  from  10  to  20  atm. 

b)  A  second  air  flowmeter  in  parallel  with  the  first  makes  pos¬ 
sible  very  accurate  air  flow  readings  in  the  low  velocity  range. 

c)  A  fourth  heater  in  series  with  the  three  already  existing 
was  added  to  raise  the  maximum  inlet  temperature  to  ^900°K. 

d)  A  sixteen-loop  gas-sampling  and  storage  valve  controlled  by 
its  own  microprocessor  has  been  connected  to  the  gas  chromato¬ 
graph.  This  valve  can  perform  automatic  sequential  GC  analysis 
of  samples  acquired  during  the  course  of  experimental  runs. 

Operation  of  the  apparatus  is  described  below.  Preheated 
air  at  a  measured  flowrate  is  supplied  to  a  690  mm  long  test 
section  with  25.4  mm  square  channel.  A  catalyst  is  placed  with 
its  downstream  end  90  mm  from  the  test  section  outlet,  and  in¬ 
sulated  from  the  wall  by  Fiberfrax  paper.  Details  of  the  cata¬ 
lyst  section  are  shown  in  Fig.  2,  A  fuel  injector  consisting  of 
five  1.6  mm  diameter  tubes,  each  containing  five  0.3  mm  diameter 
holes,  is  located  440  mm  fron  the  catalyst  inlet.  A  combination 
pitot  tube  and  thermocouple  is  mounted  200  mm  from  the  catalyst 
inlet.  In  addition  to  measuring  gas  velocity  and  temperature, 
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fig.  l  CATALYTIC  COMBUSTION  FLOW  SYSTEM 
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pitot  tube  is  used  to  extract  gas  samples  which  are  analyzed  to 
determine  equivalence  ratio.  Pressure  is  regulated  by  a  valve 
in  the  exhaust  pipe,  and  taps  placed  up  and  downstream  of  the 
catalyst  are  used  to  measure  inlet  pressure  and  pressure  drop. 

A  mass  flowmeter  (Hastings  Model  AHL-100P  with  H-3M/L-100  Trans¬ 
ducer)  measured  the  air  flowrate.  The  water  content  of  the 
inlet  air  was  measured  using  a  semiconductor  sensor  (Thunder 
Model  2000  with  BR-101B  probe)  mounted  in  the  airstream  between 
the  receiving  tank  and  the  heaters.  Inlet  conditions  for  the 
CO/air  runs  were  the  following: 

Inlet  temperature  (T^  )  =  600  ±  10  K 

Inlet  pressure  (Pi  )  =  110  ±  5  kpa  and  200  ±  7  kPa 

Inlet  velocity  (u.  )  =  10  ±  4  to  70  ±  9  m/s 

in 

CO/air  equivalence  ratio  (4>)  =  0.03  ±  0.006  to  0.32  ±  0.07 

H20  =  0.54  ±  0.007  moll 

For  the  C.jHg/Air  runs  the  inlet  conditions  were: 

Inlet  temperature  (T.  )  =  650  ±  13  K  to  800  ±  16  K 

in 

Inlet  pressure  (pin)  =  HO  ±  5  kPa 

Inlet  velocity  (uin)  =  10  ±  4  to  40  ±  7  m/s 

C-jHg/air  equivalence  ratio  (4>)  =  .19  ±  .03  to  .32  ±  .04 

H20  =  1.2  ±  .6  to  1.7  ±  .6  moll 

A  continuous  effort  has  been  made  to  minimize  crosswise 


gradients  in  temperature,  velocity,  and  fuel  concentration  of 
the  inlet  stream.  The  entire  test  section  is  insulated  so  that 
uniform  temperature  across  the  width  of  the  test  section  is  ob¬ 
tained  when  sufficient  time  (^1  hr)  has  elapsed  after  startup 
of  the  air  preheat  system.  The  uniformity  is  independent  of 
the  air  flowrate  over  the  range  of  velocities  used.  Sufficient 
fuel/air  mixing  was  obtained  only  by  placing  screens  down¬ 
stream  of  the  fuel  injector,  and  these  affect  the  velocity  pro¬ 
file.  By  trying  various  configurations,  the  velocity  unifor¬ 
mity  was  improved  while  maintaining  an  even  distribution  of 
fuel.  The  arrangement  used  in  the  present  experiments  consisted 
of  two  screens,  each  containing  four  8  mm  diameter  holes,  placed 
30  and  110  mm  downstream  of  the  fuel  injector.  Another  screen 
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perforated  with  1.6  mm  diameter  holes  was  190  mm  from  the  in¬ 
jector  and  50  mm  upstream  of  the  pitot  tube.  The  resulting  fuel 
distribution  shows  good  uniformity  (+3%)  over  the  measured  range. 
Velocity  profiles  are  less  satisfactory  (range  +6%  over  50%  of 
the  channel  width) .  Average  reference  velocities  were  determined 
using  the  CO  and  air  flowrates,  inlet  temperature  and  pressure, 
and  the  cross  section  area  of  the  catalyst.  The  catalyst  inlet 
velocity,  taking  into  account  the  fraction  of  open  monolith  area, 

is  u.  =  (1.67  +  .06)u 

in  —  rer 

Substrate  temperatures  were  measured  by  a  method  similar  to 
that  described  by  Kesselring,  Krill,  and  Kendall  (1977) .  Ni-Cr/ 
Ni-Al  thermocouples  are  fed  through  the  test  section  wall  and 
into  the  ends  of  catalyst  channels.  The  lengths  of  wire  inside 
the  catalyst  are  covered  by  mullite  insulator  and  both  ends  of 
the  channel  sealed  with  ceramic  adhesive.  The  lifetime  of  these 
thermocouples  under  test  conditions  is  short  (5-20  hr).  Ex¬ 
haust  gas  samples  were  taken  through  an  expansion  quenched, 
water  cooled,  stainless  steel  probe  mounted  in  an  elbow  down¬ 
stream  of  the  test  section.  Exhaust  gas  temperature  was  measured 
with  a  thermocouple  mounted  on  the  probe.  CO  and  CO2  were  deter¬ 
mined  by  infrared  absorption  (Horiba  Model  AIA-21) ,  C>2  by  mag¬ 
netic  susceptibility  (Scott  Model  250) ,  and  total  hydrocarbon 
(HC,  reported  here  as  C^)  by  a  flame  ionization  detector  (Scott 
Model  415) . 

The  catalyst  was  platinum  supported  on  split  cell  corrugated 

Cordierite  with  y-alumina  washcoat.  The  overall  dimensions  of 

the  monolith  were  24  x  24  x  76  mm,  the  open  area  was  64%,  the 

2 

channel  cross-section  area  was  1.9  mm  ,  and  the  platinum  loading 
3 

was  4.2  kg/m  .  The  physical  properties  of  the  catalyst  are  list¬ 
ed  in  Table  1.  The  sample  was  pretreated  by  burning  propane  for 
two  hours  with  the  maximum  substrate  temperature  at  1480K.  The 
fuels  were  carbon  monoxide,  C.P.  99.5  mole%  min,  and  natural 
propane,  96  mole%,  nominal. 
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Table  1.  Catalyst  Properties 

SUBSTRATE:  Cordierite,  American  Lava  Corp. ,  AlSiMag  795, 
split  cell 

length  0. 0760 ±  0. 0003  m 

-3 

wall  thickness  0.25x10  m 
open  area  64% 

6  —  2 

channels  per  unit  area  0.34x10  m 

—  6  2 

open  area  per  channel  1.9  x 10  m 

2  3 

ratio  surface  area  to  total  volume  2100  m  /m 

2 

ratio  surface  area  to  gas  volume  3260  m  /m~ 
channel  hydraulic  diameter  (4/3260)  0.00123  m 

3 

bulk  density  610  kg/m 

3 

solid  density  1700  kg/m 

safe  operating  temperature  1473K 

specific  heat  800  J/kg*K 

coefficient  of  thermal  expansion  (linear,  294-1033K)  3.8x10  ®  K 
thermal  conductivity  of  solid  (572K)  1.4  J/m*s*K 

approximate  channel  cross-section  is  a  trapezoid  with  base 
lengths  .001m  and  .0023m,  and  height  .0012m 

WASHCOAT:  y-alumina 

loading  115-125  kg/m^ 

surface  area  (BET)  (29.2-33.0)  x 10^  m^/m^ 

CATALYST:  platinum 

3 

loading  4 . 2  kg/m 

4  2  3 

surface  area  (CO  chemisorption)  6  x 10  m  /m 

The  y-alumina  and  platinum  were  applied  by  Matthey  Bishop,  Inc. , 
Malvern,  Pennsylvania. 
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B.  MATHEMATICAL  MODEL 


For  simplicity  only  an  individual  channel  of  the  catalytic 
monolith  was  modeled.  This  assumes  that  in  the  radial,  or  cross¬ 
wise  direction,  the  temperature  distribution  inside  the  sub¬ 
strate  is  symmetrical,  so  that  heat  transfer  between  adjacent 
channels  can  be  neglected.  In  actuality  insulation  of  the  mono¬ 
lith  and  of  the  test  section  limits  crosswise  gradients  to 
^6  K/itni,  or  ^10  K/channel.  Since  the  embedded  thermocouples 
and  the  combination  thermocouple/sample  probe  measure  quantities 
related  to  the  central  channels,  errors  due  to  nonuniform  chan¬ 
nel  performance  were  minimal.  Axisymmetry  was  assumed  in  model¬ 
ing  the  individual  channel,  whose  cross-section  is  roughly  trap¬ 
ezoidal.  The  diameter  of  the  model  channel  was  chosen  as  the 

2  x  Ayps 

hydraulic  diameter  D,  =  - rr— r — •  For  the  monolith  used  in 

J  h  perimeter 

this  study  =  1.4  mm. 

The  gas-phase  solution  inside  an  individual  monolith  chan¬ 
nel  is  based  on  the  TEACH  code  documented  by  Gosman  and  Ideriah 
(1976) .  The  general  conservation  equation  is 


TE*  +  F  [j*  (pru<l>)  +  Tr"  {prv<!)) 


_9_ 

3x 


(rr 


3|.  _  J_ 

dx  9r 


S 


<P 


where  <p  =  u,  v,  T,  Y 
property,  such  as  v, 
term,  such  as  - 
equations. 


etc.;  T  is  the  appropriate  transport 
k,  PD^K»  and  is  the  appropriate  source 
or  the  rate  of  formation  for  the  species 


Molecular  transport  properties  were  used  throughout  the 
calculations,  as  the  Reynolds  number  (based  on  inlet  flow  con¬ 
ditions)  was,  at  most,  =1900. 


a)  CO  Chemical  Kinetics  Model 


The  model  chosen  for  the  gas-phase  kinetics  of  CO  with  O, 

1  i 

is  a  simple,  one-step  overall  reaction  CO  +  j  02  ■+  COj.  The 
rate  was  taken  from  Howard,  Williams  and  Fine  (1973) : 


=  -Kq  [C0][02]°-5  [H20]0-5  exp(-E/RT) 

14  3 

with  K  =  1.3x10  cm  /mole*  sec 

E  =  126  KJ/mole  =  30  Kcal/mole 
Justification  for  this  model  lies  mainly  in  its  simplicity. 

The  model  for  the  catalytic  wall  assumes  that  the  rate  of 
CO  oxidation  at  the  wall  is  infinitely  fast,  i.e.  that  the  CO 
concentration  goes  there  to  zero.  This  was  later  confirmed 
experimentally,  and  also  is  the  only  assumption  that  could 
reasonably  reproduce  the  experimental  emission  data. 

b)  C3H8  Chemical  Kinetics  Model 

In  the  gas  phase  a  set  of  three  overall  reactions  based  on 
the  Hautman  model  (1981)  describes  pyrolysis  of  C^Hg  with  for¬ 
mation  of  an  intermediate  (C2H4) ,  oxidation  of  the  intermediate 
to  CO,  and  finally  CO  combustion  to  C02.  The  actual  reaction 
rates  in  the  present  work  are  in  the  original  form  Hautman  used 
initially  to  correlate  his  data  (Hautman,  1980) .  The  justifica¬ 
tion  for  choosing  it  is  in  its  convenience:  the  rates  stay  al¬ 
ways  finite,  while  in  the  more  recent  form  the  inverse  dependence 
of  the  C2H4  oxidation  rate  on  C^Hg  may  result  in  an  excessively 
large  rate  wherever  C-^Hg  goes  to  zero.  Because  of  the  very 
nature  of  the  numerical  solution,  this  indeed  may  happen;  to 
prevent  this  type  of  problem  the  old  form  for  the  rates  was  used. 
This  did  not  change  the  corresponding  rates  at  the  average  tem¬ 
perature  of  the  gas  inside  the  pipe. 

The  3-step  mechanism  for  propane  oxidation  in  the  gas 
phase  was  essential  for  a  good  agreement  with  the  experimental 
data.  (No  reasonable  results  could  be  obtained  by  simply  assum¬ 
ing  CgHg  to  break  down  into  CO  and  then  oxidizing  CO  to  C02) . 

At  the  platinum  surface  the  same  three  fuel  species  present  in 
the  gas  react  to  form  products.  Pyrolysis  rates  for  propane 
are  in  fact  much  slower  than  diffusional  or  convection  rates 


for  the  geometry  in  question  (Hautman  et  al,  1981).  Evidence 
for  the  formation  of  lighter  intermediates  at  the  surface  is 
scarce;  however  it  is  difficult  to  imagine  a  straight-chain 
hydrocarbon  breaking  all  its  C-H  bonds  simultaneously  to  form 
CO 2  and  H20.  In  fact,  lower  temperature  studies  show  that  pro¬ 
pane  dehydrogenates  to  propene  on  Platinum  (Biloen  et  al,  1977) . 
It  is  therefore  conceivable  that  formation  of  intermediates, 
and  possibly  radical  species,  take  place  very  close  to  the 
catalyst,  i.e.  in  a  thin  sublayer  embedded  in  the  channel  boun¬ 
dary  layer.  If  so,  only  optical  techniques  could  detect  them. 

In  order  to  simplify  the  model  both  conceptually  and  numer¬ 
ically,  the  assumption  was  made  of  one-step  reactions  at  the 
catalytic  wall.  Heterogeneous  rates  data  were  obtained  from 
Schwartz  et  al  (1971) .  As  for  CO,  inf initely-f ast  kinetics  was 
assumed.  With  these  assumptions,  the  reactions  and  their  rates 
are  shown  in  Table  2. 

To  test  the  validity  of  the  model  and  separate  gas  phase 
effects  from  substrate  effects,  the  gas  phase  equations  were 
solved  by  uncoupling  them  from  the  substrate  energy  equation. 

In  this  case  the  only  extra  unknown  variable  is  the  wall  tem¬ 
perature,  which  was  taken  from  the  experiments  and  imposed  as 
boundary  condition  at  the  wall.  This  procedure  has  as  advan¬ 
tage  the  implicit  inclusion  of  radiative  effects  in  the  model 
(provided  the  gas  is  optically  thin) ,  since  the  experimental 
wall  temperature  does  include  the  effects  of  radiation. 

C.  RESULTS 

a)  CO/Air  Mixtures 

A  total  of  78  runs  were  performed.  Analysis  of  the  Ap/p 
results  indicates  that  the  hydrulic  diameter  choice  2R=1.4  mm 
in  the  model  reproduces  the  experimental  data  (within  the  error 
limits)  at  both  p  =  110  and  p  =  200  kPa.  Trends  are  as  expected 
with  Ap/p  increasing  with  velocity  and  <f>  (i.e.  gas  temperature) 
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TABLE  II.  CHEMICAL  REACTION  RATE  DATA  USED  IN  THE 
MODELING  OF  3-STEP  KINETICS  OR  PROPANE  OXIDATION 
ON  Pt/y-Al203/C0RDIERITE  CATALYST 

REACTIONS 

GAS  WALL 


1. 

C3Hg+  j  02  -*■  |  C2H4  +  H20 

1. 

C3H8 

+ 

502  -  3C02  + 

4H20 

2. 

C2H4  +  202  ■>  2CO  +  2H20 

2. 

C2H4 

+ 

302  -*■  2C02  + 

2H20 

3. 

CO  +  02  ■+•  co2 

3. 

CO  + 

1 

7 

°2  -*■  C02 

RATES 

3  2 

GAS  (mol/cm  *s)  WALL  (mol/cm  *s) 

1.  -1.4  x  109  exp (-40608/RT)  1.  -1.51  x  103  exp (-17600/RT) 

[c3h8] 

•  (CC3H8]/[C3H8]in) -25  C02 11 * 04  2.  -80  exp(-12000/RT)  [C2H4] 

2.  -6.15  x  1013  exp (-60184/RT)  3.  -3.83  x  103  exp (-24900/RT) 

[CO] 

•  <[c2H4vrc3H8]in>-96  [o2]x-18 

3.  -2.2  x  1012  exp (-40000/RT) 

•  [C0][02r25[H20]*5 
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since  Ap/p  grows  with  both  viscosity  and  velocity  gradient  at 
the  wall  (see  Figs.  3  and  4) .  A  sample  of  wall  temperature 
profiles  is  shown  in  Figs.  5,  6  and  7.  At  the  same  equivalence 
ratio  there  is  hardly  any  effect  due  to  velocity  or  pressure, 
meaning  that  the  overall  heat  release  process  is  diffusion- 
controlled.  Flow  velocity  does  not  affect  heat  release  at  the 
wall  since  it  does  not  affect  CO  diffusional  transport  rate; 
as  for  pressure,  the  reduction  in  diffusional  transport  at 
higher  pressures  is  compensated  by  the  reduction  in  thermodif- 
fusivity  a  =  k/Cpp,  provided  of  course  that  Le  =  1. 

Emissions  in  terms  of  £  =  xcq/(xco  +  XC02^  ^*e*  terms  °f 
unburnt  carbon,  are  plotted  versus  <p  in  Figs.  8  and  9.  The 
theoretical  model  predicts  emissions  well  at  high  and  inter¬ 
mediate  velocities  while  is  less  satisfactory  when  the  velocity 
is  low.  The  main  features  of  CO  combustion  can  be  deduced  from 
Fig.  8.  At  a  constant  flow  velocity  and  -  0 ,  xcc/(xco  +  XC02^ 
is  unity.  As  <t>  is  slightly  increased,  the  catalytic  reactions, 
slow  because  of  <j>  »  0 ,  tend  to  reduce  CO,  and  emissions  drop. 

This  regime, in  which  only  part  of  the  CO  transported  at  the 
wall  by  diffusion  is  catalytically  oxidized  is  the  kinetically- 
controlled  regime,  where  finite  CO  conversion  rates  (of  the 
type  suggested  by  Voltz  et  al.  in  1973)  apply. 

Still  increasing  <)>  accelerates  the  wall  reaction  (s)  until 
they  become  as  fast,  or  faster,  than  the  rate  of  transport  to 
the  wall.  At  this  point  the  regime  becomes  diffusion-controlled 
which  is  shown  by  an  invariance  with  respect  to  <J>.  This  flat 
part  of  the  emission  curves  is  broader  at  higher  flow  velocities 
(short  residence  time)  and  becomes  very  narrow,  i.e.  barely 
noticeable,  for  long  residence  times. 

Eventually,  as  <t>  is  further  increased  and  the  wall  becomes 
hotter,  gas-phase  chemical  reactions  become  faster  and  more  and 
more  CO  is  oxidized  in  the  gas  phase,  bringing  down  emissions. 


PRESSURE  DROP  Ap/p  IN  THE  MONOLITH 


L2 


FIG.  3.  Pressure  Drop  Inside  the  Monolith  As  a  Function 
of  Inlet  Parameters  for  CO/Air  Oxidation  on  Pt. 
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Increasing  or  reducing  the  flow  velocity  alters  quantita¬ 
tively  but  not  qualitatively  this  conceptual  description.  In 
particular,  the  emission  plateau  corresponding  to  the  pure 
diffusion-controlled  regime  is  very  narrow  at  low  velocities 
and  broadens  as  the  velocity  is  increased,  until,  when  the 
velocity  is  infinitely  large,  occupies  the  entire  4>  range, 

i,e*  xcc/ ^ XC0  +  xco2 ^  =  1* 

The  effect  of  pressure  can  be  observed  in  Fig.  9.  Since 
the  diffusion  coefficients  are  inversely  proportional  to  pres¬ 
sure,  increasing  p  from  110  to  200  kPa  raises  emissions  and 
widens  the  purely  diffusion-controlled  region.  For  <p  2  .  3  the 
experimental  data  show  a  steeper  slope  for  p  =  200  kPa.  This 

is  again  the  effect  of  gas-phase  kinetics,  which  is  roughly  pro- 
2 

portional  to  p  ,  thus  favoring  emission  reductions  when  the 
pressure  is  raised. 

The  theoretical  model  underpredicts  emissions  at  larger  $ 
and  for  longer  residence  times,  i.e.  when  gas-phase  kinetics 
contributes  substantially  to  CO  conversion.  This  observation 
leads  to  suspect  the  overall,  one-step  CO  oxidation  rate  of 
Howard  et  al.  The  rate  was  proposed  for  lean  CO/air  reactions 
at  p  =  1  atm. ,  and  may  not  be  applicable  to  higher  pressures 
and  larger  <f>. 

An  important  question  for  practical  applications  of  cata¬ 
lytic  combustion  concerns  the  relative  magnitude  of  gas-phase 
versus  wall  CO  oxidation. 

Figure  10  shows  that  gas  phase  becomes  important  for  CO 
conversion  only  toward  the  outlet  of  the  channel,  where  tempera¬ 
ture  also  is  highest.  Integrating  the  local  wall  and  gas-phase 
contribution  over  the  entire  channel  gives  the  overall  effect 
of  the  catalyst,  which  is  dominant  at  high  velocities  and  lower 
4>.  When  residence  time  is  lengthened  and  <}>  is  higher  gas-phase 
CO  conversion  is  no  longer  negligible  (Fig.  11) . 
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FIG.  10A.  Gas  Phase  vs.  Surface  Contribution  to  CO/Air 
Oxidation  on  Pt-Effect  of  Velocity 
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FIG.  10B.  Gas  Phase  vs.  Surface  Contribution  to 
CO/Air  Oxidation  on  Pt-Effect  of 
and  Pressure 
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FIG.  11B.  Gas  Phase  vs.  Surface  Contribution  to  Co/Air  Oxidation  on  Pt 
Effect  of  <p  and  Pressure 
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b)  C0/Inert/02  Mixtures 

To  verify  the  effect  of  diffusion  on  CO  conversion  a  total 
of  51  runs  were  performed  with  fuel/oxidizer  combinations  in 
which  air  was  replaced  by  mixtures  of  an  inert  (N2,  C02 ,  He,  A) 
and  02 .  The  experimental  procedure  in  these  runs  consisted  in 
bringing  the  temperature  of  the  test  section  and  connected 
hardware  to  the  desired  steady-state  conditions  using  ambient 
air  from  the  compressor,  switching  to  the  C0/Inert/02  mixture 
for  the  actual  data  taking. 

The  influence  of  diffusion  can  be  estimated  in  terms  of  the 
binary  diffusion  coefficient,  which  for  a  mixture  of  gases  de¬ 
pend  on  the  molecular  weights  of  the  individual  components  in  a 
rather  complex  way.  For  simplicity,  since  the  mass  fraction  of 
the  inert  gas  dominates  over  all  others,  and  since  the  molecular 
weights  of  CO  and  02  are  similar,  the  effect  of  replacing  N2  by 
A,  He  or  C02  may  be  studied  by  looking  only  at  the  change  in  the 
binary  diffusion  coefficient  of  CO  and  the  inert  gas. 

For  gas  1  i  CO  and  gas  2  =  inert,  the  ratio  D1  2^D1  N  ^s: 

Inert  Gas  2  D,  _/D,  v, 

C02  0.86 

n2  1.00 

A  1.03 

He  2.84 

i.e.  the  mass  transfer  of  CO  in  He  should  proceed  approximately 
2.8  times  faster  than  for  CO  in  nitrogen. 

When  the  input  parameters  such  as  inlet  temperature,  velo¬ 
city  and  equivalence  ratio  are  kept  unchanged  and  only  the  inert 
gas  is  varied,  the  results  show  only  the  effect  of  diffusive 
mass  and  heat  transfer,  since  gas  and  surface  kinetics  are  not 
directly  affected  by  the  presence  of  the  inert  gas. 
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With  reference  to ' the  nomenclature  in  Fig.  2,  the  param¬ 
eters  and  conditions  kept  fixed  were: 


Inert  Gas: 

A 

He 

N2 

C°2 

Pressure  p  : 

o 

110  +  4 

110  +  4 

110  +  4 

110±4 

kPa 

Inlet  Temperature  Tq: 

600±10 

700+20 

700± 10 

690+14 

K 

Inlet  Velocity  u2 : 

11±6 

U> 

If 

13±5 

13+4 

m/s 

Stoich.  Fuel/ (02+Inert) 
ratio:  * 

0.297 

1.175 

0.274 

0.395 

CO  conversion  is  illustrated  in  Figs.  12,  13  and  14.  Fig. 

12  shows  the  dependence  on  equivalence  ratio  for  the  case  of 
C02  as  an  inert.  Comparing  the  results  with  the  ones  for  N2,  A 
and  He  in  Figs.  13  and  14,  it  is  clear  that  unburned  CO  increases, 
increasing  the  average  molecular  weight  of  the  gas  mixture. 

Figs.  15  through  18  show  the  experimental  measurements  of 
substrate  temperature  versus  axial  position  for  the  four  inert 
gases  tested  and  for  various  equivalence  ratios.  From  these  data 
is  also  clear  that  the  average  substrate  temperature  increases 
with  the  magnitude  of  the  binary  diffusion  coefficient  of  the 
(CO/inert  gas)  mixture  as  can  be  seen  from  the  following  table: 


Average  Substrate  Temperature,  K 


Inert  Gas 


co2 

N2 

He 

0.1 

850 

950 

1050 

0.2 

1000 

1200 

1400 

0.38 

1250 

1450 

no  data 

The  data  on  Argon  are  not  included  here  since  the  lower  inlete 
temperature  makes  unfair  the  comparison  with  the  other  three 
inert  gas  cases.  No  runs  were  attempted  with  He  as  the  inert 
gas  for  <J>  >  0.25,  since  the  cordierite  substrate  could  not 


/<xrn+x 
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p0  =  110  ±  4kPa 
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fig.  12  FRACTION  OF  CO  UNBURNED  VS.  EQUIVALENCE 

RATIO:  C0/02/C02 
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fig.  13.  FRACTION  OF  CO  UNBURNED  VS.  EQUIVALENCE  RATIO 
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withstand  the  corresponding  high  temperature.  These  results 
again  confirm  the  diffusion-controlled  nature  of  the  conversion 
process. 

Two  features  of  these  experimental  runs  are  immediately 
discernible  when  one  compares  the  runs  with  N2  as  inert  (Fig. 

16)  to  the  runs  with  air  (Fig.  8)  for  the  same  equivalence 
ratios  and  velocities:  1.  CO  emissions  are  about  10  to  100  times 
higher  and  2.  they  increase  with  tj>  rather  than  decreasing. 

Notice  that  emissions  are  higher  for  the  C0/N2/02  case  even 
though  the  inlet  temperature  is  100K  higher  than  for  the  CO/air 
case. 

Degradation  of  the  catalyst  was  at  first  suspected  and  then 
eliminated  since  the  runs  with  C0/N2/02  were  performed  immediate¬ 
ly  after  the  runs  with  CO/air.  The  only  difference  between  the 
two  series  of  runs  is  the  water  content  of  the  inlet  gas.  In 
the  CO/air  runs  the  water  molar  fraction  was  0.0054.  In  the 
runs  with  CO/Inert/02  the  water  content  was  determined  by  the 
H20  present  in  the  gas  bottles,  typically  ^0(1)  ppm,  and  by  the 
residual  water  left  in  the  test  rig  by  the  preheating  with  air. 

No  reliable  measurements  of  H20  were  taken  (due  to  the  hygro¬ 
meter  characteristics)  except  for  C0/N2/02,  where  the  H20  molar 
fraction  at  the  hygrometer  station  decreased  during  consecutive 
runs,  i.e.  as  $  was  increased,  from  0.002  to  0.0004.  It  is 
conceivable  that  H20  concentrations  may  have  been  even  lower  at 
locations  further  downstream  of  the  hygrometer  and  close  to  the 
test  sections.  The  gradual  "flushing  out"  of  H20  during  opera¬ 
tion  with  dry  mixtures  parallels  the  abnormal  emission  trend 
with  <p. 

That  CO  oxidation  depends  on  water  content  would  not  be 
surprising  at  all  if  took  place  in  the  gas  phase.  However,  both 
experimental  data  and  theoretical  predictions  indicate  that 
almost  all  CO  is  converted  to  C02  at  the  catalyst  surface  and 
not  in  the  gas,  except  for  low  velocity  and  high  $>,  where  the 
gase  phase  contribution  is  (at  most)  30%  of  the  total  CO  con¬ 
verted.  This,  however,  cannot  explain  the  order (s)  of  magnitude 


mi  i  iiiiMkWBtfiiigriiifi  '■irrlaiii'iniii 


increase  in  the  emissions.  The  only  possible  explanation  must 
be  that  1^0  has  a  strong  influence  on  the  catalytic  oxidation 
of  CO  on  platinum. 

This  conclusion  opens  the  door  to  some  interesting  specula¬ 
tions:  it  is  known  that  gas  phase  oxidation  of  CO  requires  OH 

radicals  (Gardiner  and  Olson,  1980) .  It  is  also  known  that 
platinum,  and  other  Group  VIII  metals,  tend  to  decompose  H2O 
and  H2O2  through  the  tendency  of  Pt  to  chemisorb  H  atoms.  It 
is  therefore  conceivable  that  at  high  temperature  the  conven¬ 
tional  CO  oxidation  mechanism  on  Pt,  proceeding  via  chemisorp¬ 
tion  of  either  or  both  CO  and  O2  (Engel  and  Ertl  ,  1979;  Camp¬ 
bell  et  al,  1980;  Strozier,  1980)  gives  way  to  a  different 
mechanism  in  which  OH  radicals  produced  at  the  platinum  surface 
by  the  H2O  present  in  the  incoming  stream  attack  the  CO  mole¬ 
cules  in  a  thin  gas  phase  layer  adjacent  to  the  catalyst,  thus 
actually  simulating  an  heterogeneous  oxidation.  It  is  also 
interesting  that  high  temperature  OH  desorption  from  Pt  has 
been  proposed  as  a  "clean"  source  of  OH  for  chemical  kinetics 
studies  (Talley  et  al,  1979;  Talley  et  al,  1981).  Low  pressure 
studies  indicate  that  OH  escapes  the  surface  only  above  a  tem¬ 
perature  ^800°K  (Tevault,  Talley  and  Lin,  1980) .  Even  though 
catalytic  processes  generating  radicals  have  been  observed  before 
and  the  term  heterohomogeneous  catalysis  has  been  coined  to 
describe  them  (Satterfield,  1980) ,  this  is  the  first  time  that 
they  are  reported  as  possibly  present  in  a  practical  combustion 
device.  Further  work  is  in  progress  to  check  quantitatively  the 
oxidation  mechanism  postulated  above. 

c)  C3Hg/Air  Mixtures 

Some  of  the  theoretical  results  are  shown  in  Figs.  19 
through  30.  Figure  19  shows  the  monotonical  decrease  in  HC 
concentration  along  the  axis  and  from  the  axis  to  the  catalytic 
wall.  The  rate  of  HC  breakdown  and  oxidation  is  a  function  of 
the  radial  mass  transfer  and  of  the  CO  oxidation  rate.  Increasing 
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the  inlet  velocity,  as  in  Fig.  20,  results  in  a  correspondent 
decrease  in  residence  time  and  a  slower  HC  disappearance. 

The  large  gradients  at  the  wall  and  throughout  most  of 
the  channel  cross-section  indicate  that  the  two-dimensional 
modeling  of  the  fuel  conversion  is  necessary  for  its  adequate 
description. 

CO  formation  and  destruction  is  shown  in  Figs.  21  and  22 
for  the  same  two  cases  of  low  and  high  inlet  velocity.  At 
ure£  =  6  m/s  the  residence  time  is  sufficiently  large  for  CO 
to  form  and  be  oxidized  mostly  within  the  channel  length.  The 
radial  CO  profiles  of  Fig.  22  indicate  an  initial  bulge,  or 
maximum,  for  r/R  'v  0.6  near  the  inlet,  where  gas  phase  produc¬ 
tion  of  CO  becomes  larger  than  its  consumption  both  in  the  gas 
phase  and  at  the  wall.  Further  inside  the  channel  propane  and 
ethene  pyrolysis  keep  producing  more  and  more  CO;  a  peak  appears 
in  fact  at  approximately  x/L  =  0.5.  Past  the  channel  half- 
length  the  high  temperature  speeds  up  the  gas  phase  rate  of  the 
CO  -*■  C02  reaction,  and  CO  decreases  rapidly.  Local  changes  in 
the  radial  profiles  slopes  reflect  the  interplay  of  the  three 
mechanisms  at  work  for  CO  -  production  and  consumption  by  gas 
phase  pyrolysis,  and  consumption  at  the  wall.  Generally  speak¬ 
ing,  near  the  wall  the  rate  of  formation  is  high  because  of  the 
high  temperature;  however  diffusion  is  also  fast  and  CO  consump¬ 
tion  is  high.  Therefore  near  the  wall  CO  concentration  is 
generally  low  and  is  accompanied  by  sharp  gradients. 

Away  from  the  wall  instead,  the  dominant  effects  are  CO 
formation  by  C^Hg  and  C2H4  pyrolysis,  and  oxidation  to  CO 2  in 
the  gas.  Local  bulges  indicate  CO  accumulation,  i.e.  that  the 
first  phenomenon  is  quantitatively  more  important  than  the 
second;  however,  contrary  to  what  happens  near  the  wall,  CO 
levels  may  be  high,  since  gas  phase  CO  oxidation  is  not  as  fast 
as  diffusion  transport  near  the  wall. 

At  the  higher  velocity  (uref  =  24  m/s,  see  Fig.  21)  there 
is  no  time  for  CO  to  be  oxidized  appreciably  within  the  channel. 
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and  axial  profiles  are  monotonic,  showing  that  a  CO  peak  would 
presumably  occur  downstream  of  the  outlet. 

Overall  comparison  between  theory  and  experiments  is  shown 
in  Figs.  23  through  30.  The  pressure  drop  between  catalyst  in¬ 
let  and  outlet  is  important  for  gas  turbine  applications.  There 
is  a  certain  amount  of  disagreement  between  predictions  and  data; 
this  is  likely  to  be  the  effect  of  having  approximated  the 
(roughly)  trapezoidal  shape  of  the  channel  cross-section  by  a 
circular  one.  The  trends,  however,  are  reproduced  (see  Fig.  23) . 
Figure  24  shows  good  agreement  in  the  UHC  emissions  comparison. 

All  the  trends  are  reproduced  and  do  not  present  particular 
problems  of  interpretation.  Less  quantitatively  satisfactory  are 
the  predictions  of  the  CO  outlet  emissions,  although  the  trends 
are  again  correctly  predicted.  The  increase  of  CO  with  (|>  and 
inlet  temperature  is  due  to  the  faster  rate  of  C-^Hg  and  C2H^ 
pyrolysis  forming  CO  more  rapidly  inside  the  monolith  and  even¬ 
tually  reducing  CO  emissions.  The  same  result  could  of  course 
be  obtained  by  increasing  the  residence  time,  i.e.  the  monolith 
length.  The  range  of  velocities  investigated  did  in  fact  reveal 
the  same  residence  time  effect:  the  highest  velocity  has  effect 
similar  to  the  lowest  <J>  and  T^n,  and  when  velocity  is  sufficient¬ 
ly  low  the  emission  decrease  after  reaching  a  maximum.  This 
maximum  therefore  occurs  when  the  cross-section  averaged  rate 
of  CO  production  is  balanced  by  gas-phase  and  surface  consump¬ 
tion  right  at  the  monolith  exit.  Since  these  rates  depend  on 
kinetics,  at  each  velocity  there  will  be  a  <f>  and  an  inlet  tem¬ 
perature  for  which  CO  emissions  will  peak.  Conversely,  for 

fixed  kinetics  (fixed  <p  and  T.  )  there  will  be  a  velocity  for 

in 

which  CO  emissions  will  peak. 

Reasonably  satisfactory  in  a  quantitative  sense  are  also 
C02  and  02  predictions  (Figs.  26  and  27) .  The  C02  trend  is  op¬ 
posite  to  the  trends  of  UHC  and  02,  since  C02  is  a  final  product 
of  HC  oxidation. 
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FIG.  27.  02  Outlet  Emissions  for  CgHg/Air  Oxidation 

on  Pt  -  Effect  of  <p ,  Inlet  Temperature  and 
Velocity 


Since  the  agreement  between  measured  and  predicted  quanti¬ 
ties  is  reasonable,  the  mathematical  model  may  be  trusted  to 
derive  some  interesting  conclusions  concerning  the  relative 
importance  of  the  many  physico-chemical  phenomena  that  take 
place  inside  a  monolith  channel. 

The  effect  of  diffusion  was  tested  by  varying  the  tempera¬ 
ture-dependent  fit  of  the  C^Hg  diffusion  coefficient.  Informa¬ 
tion  on  C^Hg  diffusion  is  sketchy  at  the  temperatures  of  interest 
in  the  present  work  (Lannus,  1970;  Barr  and  Watts,  1972;  the 
exponent  n  for  the  C^Hg  fit 
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was  chosen  on  the  basis  of  the  good  agreement  with  the  experi¬ 
mental  data  ar  I  of  the  theoretical  information  available.  In 
doing  so  the  effect  of  diffusion  became  readily  apparent:  pro¬ 
pane  emissions  changed  by  almost  one  order  of  magnitude  when 
n  was  changed  from  1.74  to  1.95.  Diffusion  therefore  has  a 
large  effect  on  propane  conversion.  However  a  case  for  the 
process  to  be  diffusion-controlled  may  not  be  made,  since  the 
propane  heterogeneous  rate  is  finite.  In  fact,  no  reasonable 
results  could  be  obtained  when  the  wall  reaction  was  assumed 
infinitely  fast  (i.e.  diffusion  was  assumed  ccntrolling) .  A 
further  look  at  Figs.  28  through  30  illustrates  the  situation. 
In  all  thiee  the  ratio 


/v  i  k  dV 

rGS  =  7  •  TT  (k  =  C3H8 '  C2H4 '  C0) 

U  “Wk  ^ 

is  plotted  as  a  function  of  the  inlet  parameters.  The  term 
upstairs  is  the  overall  rate  of  conversion  of  fuel  species  k 
due  to  homogeneous  reactions  in  the  gas,  while  the  term  down¬ 
stairs  expresses  the  corresponding  catalytic  rate  of  conversion 
at  the  wall. 
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For  propane  (Fig.  28)  rGg  increases  with  $  and  T^n,  i.e. 
with  chemical  kinetics.  This  means  that  the  homogeneous  conver¬ 
sion  becomes  more  and  more  important  than  heterogeneous  conver¬ 
sion.  The  effect  cf  4>  and  T^n  is  felt  both  by  and 
however  the  Lewis-Semenov  number  for  C^Hg  in  air,  i.e. 

Le  _  _ K _ 

Cp  p  DC gHg/Air 

varies  between  1.5  and  1.3  between  900  and  1200°K,  i.e.  in  the 
range  of  temperatures  of  interest  in  this  study.  Thus  the  ef¬ 
fect  of  raising  heterogeneous  kinetics  is  more  than  compensated 
by  the  larger  heat  feedback  from  the  wall  into  the  gas,  with  its 
corresponding  effect  on  homogeneous  conversion.  This  one  tends 
to  decrease  with  reduced  residence  times,  while  heterogeneous 
conversion  is  practically  independent  on  velocity,  since  the 
wall  temperature  depends  practically  on  $  and  T^n  only.  This 
explains  the  rGS  trend  with  increasing  reference  velocity.  Far 
more  interesting  is  the  actual  magnitude  of  rGg,  showing  that 
gas  phase  dominates  over  surface  in  converting  propane  (and 
C2H4  and  CO  as  well) .  This  result  should  not  be  interpreted  as 
meaning  that  the  catalytic  surface  has  no  effect  on  C^Hg  conver¬ 
sion.  In  fact  C,H0  pyrolysis  times  at  the  inlet  temperatures 

J  a  1 
of  this  study  are  exceedingly  long,  i.e.  of  the  order  of  10 

_3 

seconds,  while  residence  times  are  'vlO  seconds.  No  measurable 
propane  conversion  can  take  place  within  a  monolith  of  the 
length  used  without  a  catalyst.  The  effect  of  platinum  is 
therefore  to  release  heat  that  feeds  back  into  the  gas  at  a 
rate  faster  than  the  rate  at  which  propane  may  diffuse  to  the 
wall  to  react  (i.e.  Le  >  1).  Under  these  conditions,  and  for 
the  channel  size  used,  the  gas  is  where  most  of  the  fuel  is 
converted.  In  the  catalytic  oxidation  of  CO  in  air  (see  sec¬ 
tion  C.a)  the  opposite  result  holds  true;  the  two  situations 
differ  in  that  propane  reacts  at  the  wall  at  a  finite  rate 
rather  than  infinitely  fast,  its  diffusion  coefficient  is  less 
than  half  of  the  CO  coefficient  at  the  same  temperatures,  and  the 
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Lewis  number  is  larger  than  one. 

The  trends  observed  in  '  ig.  2  8  are  reversed  in  the  case  of 
the  two  other  fuel  species  (see  Figs.  29  and  30) .  Both  species 
are  formed  by  previous  reactions  and  are  not  initially  present 
in  the  channel.  Their  formation  takes  place  where  the  tempera¬ 
ture  is  higher,  i.e.  near  the  wall,  which  is  also  where  they 
can  diffuse  faster  to  >-eact  on  the  catalyst.  Therefore  an  in¬ 
crease  in  the  :inetics  parameters  and  T^n)  raises  the  temper¬ 
ature  preferentially  near  the  wall.  The  characteristic  times 
for  C2Hij  and  CO  homogeneous  breakdown  become  smaller  at  a  rate 
slower  than  diffusion  times  to  the  wall. 

The  effect  of  velocity  tends  to  oppose  <J>  and  T^fi.  The  long¬ 
er  the  residence  time  the  higher  the  amount  of  formed  species 
that  may  diffuse  to  the  wall  and  be  oxidized  there.  Therefore 
it  is  to  be  expected  that  raising  the  flow  velocity  tends  to 
increase  the  importance  of  gas  phase  conversion  of  C2H4  and  CO. 

From  this  discussion  it  is  apparent  that  no  dominant  elemen¬ 
tary  mechanism  may  be  invoked  to  characterize  propane  oxidation 
in  air  in  the  range  of  parameters  explored. 

D.  NEW  CATALYST 

The  importance  of  the  solid  phase  in  the  oxidation  of  C^Hg 
and  CO  has  pointed  out  the  necessity  of  investigating  in  a  more 
detailed  physical  way  the  nature  of  the  catalytic  material.  This 
requires  an  understanding  of  the  bulk  and  surface  properties  of 
the  catalyst  and  their  temperature  dependence  as  well  as  the 
establishment  of  criteria  that  take  into  account  chemistry  and 
physics  at  the  lattice  level  and  relate  these  phenomena  to  the 
global  catalytic  combustion  reaction. 

The  platinum-coated  cordierite  used  so  far  does  not  lend 
itself  to  an  easy  characterization,  due  to  a  lack  of  knowledge 
of  the  uniformity  and  thickness  of  the  coating,  of  the  interac¬ 
tion  between  the  y-alumina  washcoat  and  the  substrate  during 


operation,  and  of  the  aging  process  when  platinum  is  lost  due 
to  its  high  volatility.  In  addition,  the  cordierite  cannot 
sustain  temperatures  higher  than  1473°K  and  this  represents  an 
undesirable  barrier  from  the  thermodynamic  viewpoint. 

In  light  of  the  difficulties  inherent  in  the  coated  sub¬ 
strates  mentioned  above,  efforts  have  gone  into  the  design  of  a 
monolithic  single-phase  catalytic  material  with  a  low  light  off 
temperature  and  an  improved  high  temperature  stability. 

Initial  considerations  of  the  nature  of  the  catalytic  pro¬ 
cess  in  which  a  gas  phase  molecule  is  adsorbed  and  undergoes 
reaction  followed  by  desorption  of  product(s),  leaving  the  ad¬ 
sorption  site  free  for  subsequent  iterations,  indicate  that  the 
surface  sites  must  be  capable  of  (reversible)  changes  in  valence 
state.  This  is  most  readily  assured  by  making  the  material  an 
electronic  conductor.  A  substantial  body  of  research  has  been 
performed  in  the  investigation  of  high  temperature  conductive 
materials  for  use  as  MHD  electrodes,  and  the  most  promising 
candidates  for  that  application  are  perovskite  ceramics.  Perov- 
skites,  of  general  formula  ABC>3,  have  a  cubic  (or  near-cubic) 
lattice  in  which  the  B  sites  occupy  the  corners  of  the  unit 
cell,  an  A  site  is  in  the  cell  center,  and  the  oxygen  atoms  are 
at  the  edge  centers.  Typical  MHD  materials  that  have  good  ero¬ 
sion  resistance  at  high  temperatures  as  well  as  room  temperature 
resistivities  in  the  range  1-10  ftcm  are  modifications  of  the 
perovskite  compounds  LaAlO^,  MgCrO^,  etc.  It  is  possible  to 
introduce  the  necessary  properties  for  enhanced  catalytic  ac¬ 
tivity  while  still  preserving  high  temperature  stability, 
electronic  conductivity  and  low  erosion  rate  by  suitable  modi¬ 
fication  of  the  A  and  B  sub-lattice  constituents  (See  Fig.  31) . 

Introduction  of  a  strontium  dopant  on  the  A  sites  in  the 
material  La  g^Sr  3gCr03  produces  a  p-type  material  having  de¬ 
creased  electrical  resistivity  as  well  as  increased  volatility 
(Meadowcraft ,  1969)  (See  Fig.  32) .  To  reduce  the  volatility, 
aluminum  can  be  introduced  into  the  B  sites,  as  in  LaAl  ^Cr  503, 
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Fig.  31.  Perovskite  Crystal  Structure  (from 
Voorhoeve  et  al.f  1976) 
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which  is  extremely  stable  at  high  temperatures.  Mass  loss  mea¬ 
surements  on  this  material  show  that  the  mass  loss  rate  declines 
with  increasing  A1  content  up  to  atomic  fraction  0.5  of  the  B 
sites,  after  which  the  loss  rate  remains  essentially  constant; 
the  aluminum  serves  to  stabilize  the  chromium  in  the  lattice, 
which  otherwise  volatilizes  as  chromium  oxides  (Anderson,  1978) 
(Fig.  33) . 

Catalytic  behavior  of  mixed  perovskite  materials  shows 
moderate  activity  for  CO  oxidation  but  high  light-off  tempera¬ 
tures  and  susceptibility  to  poisoning  by  S02>  Addition  of 
platinum  in  quantities  >1600  ppm  to  the  materials  La  ?Pb  3Mn03 
and  La  ^Sr  3MnC>3  resulted  in  a  lowered  light-off  temperature 
and  performance  comparable  to  commercial  supported  Pt  catalysts 
(Voorhoeve,  1976) . 

Motivated  by  the  general  influences  outlined  above  of 
material  composition  on  performance  characteristics,  we  have 
chosen  the  material  La [Mg  Pt  Cr  ...Al  ._]0o,  with  .001  <  x  < 
.01,  as  our  monolithic  combustion  catalyst.  We  estimate  this 
material  will  exhibit  catalytic  activity  comparable  to  that  of 
platinum  coated  substrates,  remain  as  a  single  phase  over  a  wide 
temperature  range,  and  have  suitable  low  light-off  temperatures. 
The  platinum  is  chemically  bound  into  the  material  rather  than 
present  as  islands  of  free  metal,  and  can  therefore  be  expected 

to  have  low  volatility.  The  magnesium  will  exist  in  the  lat- 

.  2+  .  4+ 

tice  as  Mg  which  will  serve  to  charge-compensate  the  Pt  and 

make  its  incorporation  into  the  lattice  more  probable. 

A  sample  of  this  material,  with  x  =  .01,  has  been  prepared 
for  us  (by  Prof.  H.U.  Anderson,  University  of  Missouri/Rolla) , 
and  an  x-ray  diffraction  analysis  (by  Dr.  G.J.B.  Williams  and 
Mr.  J.  Hurst,  Brookhaven  National  Laboratory)  of  the  powder 
indicates  that  the  material  is  approximately  cubic  and  consists 
of  a  uniform  single  phase,  with  lattice  parameter  3.855  8,  in 
the  range  expected  for  perovskites.  Atomic  adsorption  measure¬ 
ments  performed  at  Princeton  have  indicated  the  presence  of  the 


La,.,  Sr.CrO,.^ 
25°C 


FIG.  32 
from 

Meadowcraft , 
1969) 


IT2 

•?>  E 


§  V 


0  0.1  0.2  0.3 

Sr  CONTENT.  X 

Electrical  Resistivity  of  La.  Sr  CrOj  as  function  of  Sr 
content  * 


120 

E 

o> 

3 

80 

FIG.  33 

V) 

1/) 

from 

o 

Anderson, 

mJ  , 

H- 

40 

1978) 

x 

O 

UJ 

5 

Lo  j  (Mg^  Cf,.,  Al,)jj  Ojjj 
ATM:  P.  «  10*  atm 
P. .  =  10  atm 


!740eC 


H 


0.4  0.6 

AL  CONTENT,  X 


Height  loss  (Mg/c«*/h)  of  La.t(Mg.o*Cr,  Al) .«s02.(>  as 
function  of  Al  content:  Temperature  “  1740*0;  Ata:  Flowing 
gas  mixture  (0.001  ata  0i+0.1  ata  HjO+N2);  Flow  rate  “  1 
linear  ca/sec. 


expected  amount  of  platinum;  however,  since  this  quantity  is 
small  it  is  still  possible  that  it  is  present  in  a  second  phase 
which  is  not  detectable  by  the  x-ray  measurements. 

Our  design  for  the  catalyst  unit  geometrical  configuration 
is  expected  to  prove  practicable  for  actual  combustion  applica¬ 
tions,  relatively  easy  to  fabricate  and  have  well-defined  flow 
channels  amenable  to  exact  flow  modeling  calculations.  The  lat¬ 
ter  criterion  is  essential  for  model  development  but  will  not 
be  a  constraint  on  actual  combustor  design.  The  catalyst  con¬ 
figuration  will  consist  of  a  stack  of  flat  plates  separated  by 
spacers  at  their  edges,  creating  high  aspect  ratio  channels  be¬ 
tween  the  plates.  The  plate  dimensions  will  be  approximately 
2"x5"x.050",  and  the  spacing  will  be  variable  from  .050"  to 
^0.2" .  These  dimensions  have  been  chosen  to  maximize  resistance 
to  thermal  shock,  according  to  the  criteria  developed  by 
Pfefferle  (1980). 

Knowledge  of  the  physics  of  the  new  catalyst  is  essential 
for  more  realistic  understanding  and  modeling  of  the  heterogen¬ 
eous  rates.  The  steps  to  go  from  reactants  to  products  in  a 
catalytic  process  can  be  summarized  as: 

gas-phase  diffusion 
pore  diffusion 

chemisorption  on  active  sites 

chemical  reaction (s) 
pore  diffusion 
products  pore  diffusion 
products  gas  phase  diffusion 

There  is  evidence  that  at  the  temperatures  typical  of  catalytic 
combustion  the  surface  of  the  catalyst  is  sintered  to  the  point 
that  the  surface  area  available  and  geometric  area  practically 
coincide:  for  monolithic  single-phase  materials  pore  diffusion 
at  high  temperature  may  be  also  negligible;  however,  many 
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reactions  are  possible  depending  upon  the  chemisorbed  species. 
For  instance,  in  the  C0/0 2  case,  the  possible  states  of  the; 
reactants  are:  (subscript  a  =  adsorbed,  g  =  gas  phase) 
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COa 

and 

°2a 

2) 
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and 

CO 
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A  Langmuir-Hinshelwood  rate  may  be  derived  for  each  of  these 

combinations.  For  instance,  for  3)  the  rate  of  conversion  of  CO, 

r  ,  is  (Shishtf,  1972) 
c 


pco2(p”  } 


1/2 


2  (  C02  °?  \ 

V*  ]KCOK02  ho\  -  K  1*1  } 


r_  = 


(l  +  Kcopco  *  Ko/o2> 


where  KCQ,  K0  »  Kg  are  the  equilibrium  constants  for  adsorption 
of  CO,  Oj  and^the  forward  rate  of  formation  of  active  sites. 

It  is  clear  that  only  through  a  better  knowledge  of  the 
surface  can  better  rates  be  formulated.  It  is  also  extremely 
interesting  to  find  which  of  the  steps  listed  before  is  con¬ 
trolling.  Given  the  single-phase  structure  of  the  proposed  new 
catalyst,  answers  to  the  questions  posed  by  the  considerations 
above  should  be  easier  and  far  more  reliable  than  for  coated 
catalysts. 

The  actual  forming  of  the  flat  plates  and  spacers  made  of 
the  new  catalyst  was  contracted  with  Trans-Tech,  Inc.,  Gaithersburg, 
Md.,  a  ceramic  manufacturer.  Perovskite  powder  to  the  specified 
stoichiometry  was  formulated  and  fabrication  of  the  flat  plates 
started  approximately  in  December  1980.  Several  problems  were 


encountered  in  the  process:  the  most  important  were  an  unexpec¬ 
ted  high  sensitivity  to  the  temperature  and  atmosphere  composi¬ 
tion  of  the  furnace  during  the  firing  process,  and  severe  war- 
page.  The  first  problem  made  impossible  for  a  long  time  the 
attainment  of  the  desired  density  (and  consequent  mechanical 
properties) .  The  second  made  very  difficult  or  impossible  the 
final  grinding  to  the  required  geometrical  dimensions.  Solu¬ 
tions  to  these  problems  have  currently  evolved  into  a  strategy 
of  firing  in  forming  gas  at  the  elevated  temperatures  necessary 
for  perovskite  sintering  (up  to  1700  C) ,  followed  by  reoxida¬ 
tion.  Samples  prepared  in  this  manner  are  fully  dense  and  ap¬ 
pear  to  have  good  mechanical  and  thermal  stability.  The  warpage 
has  been  attributed  to  a  possible  contamination  by  silica  present 
in  the  Trans-Tech  kilns  during  firing,  and  steps  have  been  taken 
to  control  this.  At  the  end  of  the  contract  for  which  this  re¬ 
port  is  being  written,  test  plates  were  being  evaluated  for 
platinum  content  and  qualitative  light-off  behavior. 
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